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Abstract 

We study FRW cosmology for a double scalar-tensor theory of gravity where two scalar fields are 
nonminimally coupled to the geometry. In a framework to study stability and attractor solutions 
of the model in the phase space, we constraint the model parameters with the observational data. 
For an accelerating universe, the model behaves like quintom dark energy models and predicts a 
1 ^ I transition from quintessence era to phantom era. 
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1. INTRODUCTION 



The acceleration expansion of the universe is supported by observations of high redshift 



type la supernovae, the surveys of clusters of galaxies 
(SDSS) jsl and Chandra X-ray observatory 
(CMB) anisotropies observation 



2||3|j^, Sloan digital sky survey 



6[ . In addition, Cosmic Microwave Background 
7| exhibit that the universe flatness js]. The observations 
determines, with high precisions, the basic cosmological parameters and strongly indicates 
that the universe is dominated by a smoothly and slowly varying dark energy (DE). A 
dynamical equation of state ( EoS) parameter can be considered as an effective parameter 

9| 10|. The parameter is connected directly to 



to explain the acceleration of the unvierse 
the universe energy density and indicate the expansion of the universe 



9| 10[ . In particular. 



a proposal to explain the recent observations is the quintom dark ener gy, constructe d by 
two scalar fields, and its EoS parameter crosses the phantom divide line [ni Q Q Q H- 
Motivated from string theories, the scalar-tensor models provide the simplest model- 
independent description of unification theories which predict couplings between scalars and 
curvature. They have assumed a prominent role since any unification scheme, such as 
supergravity, in the weak energy limit, or cosmological models of inflation such as chaotic 
inflation, seem to be supported by them 



the current acceleration of the universe 



16 1. In addition, they have been employed to study 
17l|l8l|l9||2o||2l||23|. 



In this paper, we study the detailed dynamics of the double scalar-tensor cosmological 
models. Since the major difficulty in cosmological models is the nonlinearity of the field 
equations and thus limitation in obtaining the exact solutions, we investigate the asymptotic 
behaviour of the model, which provides the relevant features to be compared with the 
current physical data available for the universe. In this context, the perturbation methods, 
especially linear stability analysis which have been used to study the qualitative analysis of 
the equations and of the long term behaviour of the solutions are being proposed in this 

nnnnnnn 

work [2J] [25[ [26[ [27| [28[ [29[ [30| . Sec. two is devoted to a detailed formulation of the model. 
In Sec. three, we obtain the autonomous equations of the model and the late time attractor 
solutions by using the phase plane analysis in addition to best fitting the model parameters. 
We also examine the behavior of the EoS parameter of the model which predict a transition 
from quintessence era to phantom era in future. In Sec. four, we present summary and 
remarks. 
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2. THE MODEL 



A general action in four dimensions, where gravity is nonminimally coupled to two scalar 



16| 



fields with no interaction between fields, is given by 

S = J [F{^)R + G{i^)R - \p,^^^^ + V{^) (1) 

where the functions F{ip), V{lp), G^i/j), and W^ip) are not specified. By using the trans- 
formations (j) = F[(p) and uj{(j)) = 2dF{^)/dip ' Brans-Dicke action can be simply recov- 
ered. Furthermore, in our units, the standard Newton coupling is regained in the limit 
F{ip) + 0(1/)) — )■ — By varying the action with respect to the metric g^^, the field 
equations can be derived: 

[F{^) + GimR,u - \g,uR) = T(f + (2) 

where the effective stress-energy tensors for the scalar fields Lp and "0 are: 

- g^uUF{^) + F(v^);^,, (3) 

- g,,UG{i^) + G{i^).,^,. (4) 

The variations with respect to if and if) give the klein-Gordan equations 

^AF^ dV ^ 

, ^,dG, dW 

Let us now take into account a FRW metric of the form 

ds^ = de-a^[^^^+r^dn\ 

The field equations then become 

^, , a ,d,n k'^ . ^d,,dF -dG, 
{F + G 2- + i-y + -) = -2-{ip— + ij—) 
a a a a d(p dip 

dif"^ dip dip"^ dtp 

- lil^' + i^'-v-w), (7) 



,a,n a, dF :dG. 
6{F + G = -6- vi— + 

a a dip dip 

+ 6- V + G) + hip' + ij') + V + W, (8) 

a ,d ,a,n k'^.dF dV , , 

(^ + 3- = 6(- + (-f + - )— -— , 9 

a a a a dtp dip 

h ■ /d ,d.n k'^^^dG dW ,^ , 

^ + 3-^ = 6 - + + - -7- - -7-, 10 

a a a a dtp dip 

where equation ([8]) is the energy constraint corresponding to the (0,0) component of the 
Einstein equation. 



3. STABILITY ANALYSIS AND COSMOLOGICAL TEST 



In this section, in a flat FRW universe we study the structure of the dynamical system 
via phase plane analysis, by introducing the following dimensionless variables, 

{F + G)H' 12{F + G)H^' 6{F + G)m ^ ' 

and also take A = -jy^^^. Using equations (!7|)-( !T0|) . the evolution equations of these 
variables become, 

~ X ~ 6 ^ X ^ ~ ^ ^ 6" 

- ^Z(l-A)-2yX-^, (12) 

AY AX lAY 2ZY 

Y' = 2 + + Z + — 

X 3 3 X 3 

XY . ZX 2ZY , , 

+ ^-4^2 + ^ ^ (13) 



, -hXZ AZ AZ 2Z^(1-X) 2Z^ , , 

Z' = ^ ^ - AZY + AZ, (14) 

3X3 X 3 ' ^ ^ 

where prime " ' "means derivative with respect to X = ln{a). Also, the Friedmann constraint 
equation ([8]) becomes 

X + Y + Z = l. (15) 
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In term of the new dynamical variable we also have, 

H 2 2 Z{l-\) X Z 

in which can be used to evaluate the effective EoS and deceleration parameters in terms of 
the new dynamical variables. Using the constraint (fl5|) . the equations f|T2|) - f|T^ now reduce 
to the following two equations: 

2 7 (1-X-F)(1-A) bX 



X' 



(17) 



X 3 X 2 

4X^ UY 

-^ + ^-{l-X-Y){l-X)-2YX 
{1-X-Y)X 



AY AX lAY 
Y' = 2-^ + ^-^ + (l-X-y)(l-A) 

2{1 -X-Y)Y X^ 
X 

+ iy2 il-X-Y)X 2il-X-Y)Y 

3 3 3 

- A(i-x-y). (18) 

It is more appropriate to solve the dynamical equations f ll7p - fll8p than fll2l) - (fT^ . In stability 
analysis, by simultaneously solving X' = and Y' = 0, we find the critical points (fixed 
points) and study the stability of these points. In the context of autonomous dynamical 
systems, these points are exact constant solutions and are often the extreme points of the 
orbits and therefore describe the asymptotic behavior of the system. To the first orders in 
the perturbations, by substituting linear perturbations X' — > X' + SX', Y' — )• Y' + 6Y' about 
the critical points into the two independent equations ffT7|) and ffTS]) . we yield two eigenvalues 
Aj(z = 1,2). For a stable point, we require the real part of all the eigenvalues to be negative. 
We find eigenvalues for two critical points. The expression for both critical points and 
eigenvalues in our model depend on the stability parameter A, and are highly nonlinear, 
long and cumbersome. Therefore, one can not easily reveal the nature of critical points. 
To circumvent this issue in addition to physically motivate the subject, we simultaneously 
solve the equations by best fitting the stability parameter and initial conditions with the 
observational data using the method. Next, we solve the equations by best fitting the 
model with the observational data for distance modulus. 
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3.1. Observational constraints 



The distance modulus is defined as, fi{z) = 25+5 logj^Q diiz) where the luminosity distance 
quantity, (Il^z) is 



By solving the equations, we find H{z) which can be used to evaluate fi{z). To best fit the 
model for the parameter A and the initial conditions Y{0), X{0), H{0) with the most recent 
observational data, the Type la supernovea (SNe la), we employe the distance analysis. 
We constrain the parameters including the initial conditions by minimizing the function 
given as 



where the sum is over the SNe la sample. In the above relation, and /i° are the 
distance modulus parameters obtained from model and observation, respectively. Also, the 
estimated error of the /if" is a. For our model the best fit values for the parameter and 
initial conditions occur at A = 2.02, X{0) = -0.3, Y{0) = -0.65 and H{0) = 0.91 with 
Xmin = 552.7337761. Fig. 1) shows the 1-dim and 2-dim likelihood for the parameters A 
and the Hubble parameter H{0) at the 68.3%, 95.4% and 99.7% confidence levels. 

In Fig. 2, the distance modulus, ij,{z), in our model is compared with the observational 
data for the obtained parameters and initial conditions using method. As can be seen the 
best fitted parameters and initial conditions are in good agreement with the observational 
data. In the following we investigate the stability of the model with respect to the best 
fitted model parameter. 

3.2. Stability and phase space 

Solving the stability equations for the best fitted model parameter A we find two fixed 
points with the stability properties given in tables I. 

From the above table we find that both critical points are stable. In Fig. 3, the trajectories 
entering the stable critical points FPl and FP2 in the phase plane are illustrated. For the 




(19) 




x'5ive(A,x(o),y(o),ff(o)) 



(20) 
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Fig. 1: The graph of 1-dim and 2-dhxL hkehhood distribution for parameters A and H(0) 




Bestfitof SNIa 



Fig. 2: The best fitted distance modulus fi{z) plotted as function of redshift 



given initial conditions on X, Y between -1 and 1, the trajectories (green curves) approaching 
the stable critical point FPl in the phase plane are shown. For the initial conditions on X, 
Y less than -1 or greater than 1, then the trajectory shown entering the stable critical point 
FP2. The best fitted trajectory ( red curve) with the properties given in the previous section 
enters the stable critical point FPl. This trajectory both fits the model parameter A and 
the initial conditions for X, Y and H. 
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TABLE I: Best fitted critical points 
Critical point {X, Y) stability 



FPl 



(1.72, 0.55) stable 



FP2 (-0.5,1.26) stable 



Y 0- 




/ / '/ 
^y / / / / 
y / / / iH 



X 

Fig. 3: The attractor property of the dynamical system in the phase plane. The red trajectory 
approaching the FPl corresponds to the best fitted stability parameter and initial conditions. 

3.3. Phantom crossing 



The EoS parameter is one of the most important cosmic parameters that shows the 
behavior of the universe and its dynamics. In this work, with stabihty analysis in addition 
to constraining the model with the observational data we obtain a better understanding of 
the cosmological paraments such as EoS parameter. The effective EoS parameter is given 
by uJeff = ^1 ~ 1;^ where in terms of best fitted model parameters is given in equation 
(fT6|) . From Fig. 4, the model does not predict a matter dominated epoch in the higher 
redshift. However, it exhibits phantom crossing behavior in the future (Fig. 4). The current 
best fitted effective EoS parameter is Wg// — —0.45. Fig. 4 also illustrate that the phantom 
crossing occurs in future at z ^ —0.4 and approaches —1.1 at z = —1 when a big rip occurs. 
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Fig. 4: The best fitted effective EoS parameter, We//, as function of redshift. 
4. SUMMARY AND REMARKS 



This paper is designed to study dynamics of the double scalar-tensor gravity cosmology 
by using the stability analysis and the 2-dimensional phase space of the theory. Different 
from the work in 2^, in a remarkable approach in stability analysis, we solve the system of 
autonomous differential equations by constraining the model parameters and also the initial 
conditions with the observational data for distance modulus. As a result, the critical points 
which exhibit states of the universe are observationally verified. From analysis, the two 
stable critical points in the model are depicted in Figs. 3. With the constrained parameters, 
the best fitted trajectory approaches only the first stable point. 

We then test the model against observational data by calculating the best fitted effective 
EoS parameter, Ue/f for the model in terms of the stability parameters. The result shows 
that our model exhibits a quintom behavior; an accelerating universe with a transition from 
Ueff > — 1 (quintessence era) to u^jf < — 1 (phantom era)in the future. 
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